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ABSTRACT 


Orifice and louver sound power spectra are investigated, using an 
intensity probe, at various wind speeds in a low noise, semi-anechoic, 
subsonic wind tunnel for free stream velocities below 50 meters per 
second. The radiated noise is created by turbulent flow over various 
orifice and louver geometries which are flushed mounted into the wall 
of a long duct. Five orifice samples of rectangular shape and various 
transverse dimensions as well as four louver samples with multiple 
rectangular and circular orifices are tested. Also investigated is 
the effect of the leading and trailing edge angle on the radiated 
sound power. The scaling laws of the excitation frequencies and the 
speed/power laws are presented for ratios of the boundary layer thick- 
ness to the transverse orifice dimension from 1.01 to 4.29. A 
detailed theoretical model is developed for rectangular shaped aper- 
ture orifices and louvers based on the work by Ffowcs Williams [9], 
Nelson [11], and Corcos [12]. The model describes the predominant 
variables which effect the radiated flow noise. To ascertain the val- 
idity of the model, it is used to collapse the experimental power 
spectra for two orifice and one louver geometry. The model showed 
excellent agreement with the actual sound power measurements. Direc- 
tivity plots are also presented to further describe the orifice/louver 
source. 


Thesis Supervisor: Professor Patrick Leehey 
Professor of Mechanical and Ocean Engineering 
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1. INTRODUCTION 


1.1 Background 


Turbulent flow over an orifice or louvered openings can generate 
high levels of broadband noise. This phenomenon has been known for 
some time [1] and early investigators like Dunham [2], Rossiter [3] 
and East [4] attempted to identify and control the high vibrational 
levels and possibly the radiated sound levels caused by the interac- 
tion of the flow over these openings. Recent work by Elder, Farabee 
and DeMetz [5] concentrated on the area of cavity resonance where the 
opening provides the feedback path between the closed acoustic cavity 
and the source, the turbulent flow. When the typical dimension of the 
cavity is large compared to the acoustic wavelength, the source can 
excite various room modes of the cavity. If the turbulent vortex 
shedding frequency coincides with one of these room modes, then a 
strong feedback will occur, resulting in an acoustic resonance of that 


mode. 


If the cavity's volume is further increased, multiple room modes 
will be excited and the strong feedback mechanism will be lost. This 
results in a much broader acoustic response and at lower amplitude 
levels than the resonant case which is typically characterized by 


large amplitudes at very discrete tonal frequencíes. In practice, the 
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detection of resonating cavities is often easy due to their strong 
discrete tones. However, once these resonances have been corrected, 
or presumably prevented in the design, identifying the cause of the 
remaining flow noise becomes much more difficult. Possible contribu- 
tors to the broadband flow noise are: (1) the turbulent boundary layer 
direct transfer of energy to the structure of the enclosed acoustic 
space, (2) the power flow through the opening (open window effect) 
caused by the acoustic field created by the upstream and downstream 
turbulent sources in the boundary layer, (3) the power flow into the 
acoustic space caused by the source created in the opening due to the 
driving aerodynamic field and the shear-layer edge interaction with 
the downstream edge of the orifice and (4) the vibrational excitation 
of the internal structure caused by the orifice source. The amount of 
noise which the latter contributes depends on the details of the 
internal structure and the orifice/louver attachment design.  There- 
fore to predict the contribution from internal structure vibration, 


one must first determine the orifice source. 


1.2 Content 


The contribution of the first three possible contributors to the 
total power radiated into a closed acoustic space is analyzed based on 
the results from several experiments. It is shown that indeed the 


orifice source provides both the spectral character and the predom- 
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inant amount of the total sound power. By analyzing and comparing the: 
results with a derived equivalent source model, the nature of the 
source and the parameters which affect its strength and spectral char- 
acter are identified. The analysis comes directly from experimental 
measurements taken in a subsonic wind tunnel at the M.I.T. Acoustic 


and Vibration Laboratory. 
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2. SOUND INTENSITY MEASUREMENTS 


2.1 Why Measure Sound Intensity? 


To describe an acoustic source, one must measure a quantity which 
is not dependent on the environment or the surroundings. The quantity 
most often identified to describe a source is its power. Knowing the 
source's power, its interaction in any known environment or any known 


surrounding can be determined. 


Until recently, most all experimental measurements to describe an 
acoustic source were accomplished by measuring sound pressure. This 
quantity, however, was the result of a source’s interaction with a 
specific environment and utilizing the experimental results in other 
applications was often difficult. One technique that was developed to 
overcome the difficulty in relating one environment to other was to 
isolate the source of interest and to simulate a free field environ- 
ment in an anechoíc chamber. Here the radiated sound pressure in air 
was measured, then related directly to the sound intensity. This 
technique has as its biggest limitation, the availability of an effec- 
tive anechoic chamber, particularly for use in low frequency analysis 


at the low frequencies. 


Another technique is to measure the intensity directly utilizing a 





er: 


sound intensity probe. The theory of intensity measurements comes 
directly from two fundamental properties of acoustics. Those are the 
conservation of momentum and the definition of intensity. Conserva- 


tion on momentum states that: 


pov,/dt = - др/дт (2.1) 


which can be solved for the radial component of the particle velo- 


city, у, 


vr = - (1/p) S ap/ar at (2,12) 


A finite difference approximation is used to calculate the pres- 


sure gradient using two identical microphones by 


Op/dr = lim (Ap/Ar) = (pa-pp)/Ar (2.2) 
Ar-0 
where A and B correspond to channels A and B for the respective 
microphone inputs and Ar is the finite microphone spacing. The 
instantaneous pressure is taken as the mean value of the two pressure 


signals, 


P = (Patpp)/2. (2.2a) 


The finite difference approximation is valid provided the separa- 
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tion distance, Ar, is small compared to the acoustic wavelength. 
Thus, a finite difference approximation error becomes significant at 
high frequencies. For a plane sinusoidal wave, propagating along the 
axis joining the two microphones, the finite difference approximation 
assumes the free field phase between the two microphone positions to 
be equal to kAr, whereas the actual measured free field phase 15 
Sin(kAr). Therefore, the measured intensity, I', is related to the 


actual intensity, I, by 


II -oSin(kAr)/kAr (23) 


and the resulting measured intensity underestimates the actual 


intensity by Le at high frequency by 


La, - 10 log(Sin(kar)/kAr) . (2.3a) 


All the intensity measurements in this thesis were obtained with a 
spacing of Ar = 50mm. Thus the finite difference error for this 
spacing, at the upper analysis bandwidth frequency of f = 1600 Hz, was 


Lo = - 1.7 dB. 


From the definition of intensity and by substituting equations 2.2 


and 2.2a, the measured sound intensity becomes 





I, = Pvr = - [(patpp)/2rAr)S(pa-pp)dt. (2:035) 
By measuring the intensity over a control surface which surrounds 


the source, the power can be calculated from: 


P= Jf Тү дА (2.5) 
А 


The advantages of measuring intensity directly are: (1) the 
effects from diffuse noise contaminating the experimental results are 
eliminated since the integral of this noise intensity around any 
closed control surface is zero, (2) the sound power calculated 
describes the acoustic source in a form which is readily usable to 
design engineers, (3) measurements can be made in the source's near 
field and (4) there are no restrictions on the shape and size of the 
control surface. A simple yet effective analogy to describe the 
rationale in measuring sound intensity, thus power directly follows: 

In attempting to design a heating system, an engineer 
solicits specifications from two companies. One com- 
pany describes its heating system by stating the 
electrical power and the thermal efficiency of the 
unit. The other company describes the temperature at 
various locations inside an experimental enclosure at 
Steady state conditions. 

It is obvious in this example that describing the temperature 
field within one environment is practically useless to the design 
engineer whereas knowing the power and thermal efficiency, the 


designer can determine the interaction of this source with his partic- 


ular environment. 
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In the experiments which follow, every attempt was made to uti- 
lize the advantages of measuring intensity directly and simulating a 
free field environment by taking the measurements in a semi-anechoic 
chamber, thus minimizing any enclosure interactions which could affect 


the radiated sound power. 


2.2 Ideal Versus Practical Sound Intensity Analyzer 


2.2.1 Active and reactive sound fields 


Sound fields exist in two parts: an active part which transports 
sound energy by virtue of the fact that the sound pressure and the 
particle velocity are in phase and a reactive part which stores sound 
energy and in which the sound pressure and particle velocity are in 
quadrature. A sound field is said to be more reactive the more energy 
is stored relative to that transported. Sound field reactivity is 
defined as the difference between the meaSured sound pressure level, 


Lp: and the maximum intensity level L;. 
2.2.2 Ideal sound intensity analyzer 
An ideal sound intensity analyzer would measure only the active 


part of the total sound field and would indicate sound intensity val- 


ues approaching 0 watts/m? in a highly reactive (highly diffuse) 





sound field. 


2.2.3 Practical sound intensity analyzer 


A practical sound intensity analyzer can never display such a low 
intensity level due to the small but ever present phase-mismatch 
between the two microphone channels. It is this phase mismatch which 
accounts for the lower limit in the dynamic range of a practical sound 


intensity analyzer. 


2.3 Determining the Limits of Error for the Measured Intensity 


Errors are introduced in any practical intensity measurement by: 
(1) the influence of the phase mis-match in the probe, $4 and 
the analyzer, Ha, 
(2) the influence of the phase of the sound field of known reacti- 
vity, ég, and 


(3) using a particular spacer between the two microphones. 


The limits of these errors can be determined using a phase- 
reactivity nomogram [see reference 6]. This nomogram graphically 
determines, $g for a given measured reactivity (Li-Lp) and a 


particular spacer, Ar, according to the relationship: 





DS (L5*0.16) - lOLog[1000/Ar x C/f x $£g/360] (ОХО) 


The phase of the sound field, ff, together with the phase of the probe 
and the analyzer, $f - $a - Óm and de + $a + pm» set the limits of 
error in the measured phase at a particular frequency. From these 
phase limits, the nomogram provides the limits of error in dB for the 


measured intensity. 


2.4 Sound Intensity Probe 


The specific details on the operation of the intensity probe used 
in this experiment are contained in reference (7). Briefly, the 
intensity probe consists to two B&K phase matched, 1/2" condenser 
microphones, type 4177/4165, separated by a spacer 50 mm in length and 
on a common axis. One microphone is a reference while the other 


microphone provides the pressure gradient measurement. 


2.5 Intensity Measurements 


2.5.1 Calculating intensity 


Sound intensity measurements were analyzed using a B&K 2032 Spec- 


trum Analyzer. The sound intensity was calculated according to the 


relationship: 





I(k) - - (1/oArg) Im([Gap(k)] (2.6) 


where Ar is the microphone spacing in meters, A corresponds to 
Channel A which is the reference microphone input, B corresponds to 
Channel B which is the second microphone input, k is the frequency 
index and p is the density of air, calculated according to the rela- 


tionshíp: 


p = 0.348534 p/(273+T) (2.7) 


where p is the pressure in mbar and T is the temperature in degrees 


centigrade. 


The cross spectrum, Gap(k), is the complex ensemble average of the 
complex products of the complex conjugated one-sided instantaneous 


spectrum G,(k) and the one-sided instantaneous spectrum Gp(k). 


САВ(К) = Ga(k)*Gp(k) (2.8) 
The instantaneous one-sided spectrum for channel A or B, Ga(k), 


Gp(k) was calculated according to: 


SACK) for k = 0 
Ga(k) = 2 SACK) for 1 < k < N/2 - 1 (2-9) 


0 for N/2 <k < N-l 





where N is the number of samples per ensemble. For this experiment, 
N was equal to 1024. $д(К) 15 the two-sided instantaneous spectrum 
for channel A which is determined from the recorded time record, a(n) 


according to: 


Sa(k) -F[w(n)-a(n)] (2.10) 


where 7 is the forward discrete fourier transform and w(n) is the 
time record weighting function. In this experiment, a Hanning window 


weighting function was utilized. 


The forward discrete fourier transform, 3 , is defined as: 


N-1 
Sa(k) = 1/N E a(n)+exp(-j 2rkn/N) (2.11) 
n=0 


where n equals the time index times the sampling interval, nAt 
апа k equals the frequency index times the frequency resolution, kAf. 
The frequency resolution for this experiment was 2Hz for an analysis 


bandwidth of 1600 Hz. 


2.5.2 Measurement techniques 
Two techniques can be used to measure sound intensity using an 


intensity probe. They are: 





(1) continuous spatial averaging 


(2) discrete point averaging 


Continuous spatial averaging consists of sweeping the intensity 
probe at a controlled rate over each surface while continuously samp- 
ling and averaging data. This technique produces the most accurate 
spatially averaged results but suffers slightly in repeatability. 
Several repeated runs, often four or more, were required to obtain the 
required repeatability of total sound power (within 0.5 dB on at least 
three runs). Figures 2.1, 2.2 and 2.3 are the results of three such 
repeatability experiments with 0.5 dB or less variation in total sound 
power. Shown is the sound power spectral level from a rectangular 


orifice, 4 1/2" by 1 2/3", and an analysis bandwidth of 2 Hz. 
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Figure 2.1: Louver L1 sound power spectra repeatability. Shown are four succes- 
sive spectral measurements of sound power. 
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Figure 2.2: Orifice 01 sound power spectra repeatability. 
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Figure 2.3 Orifice 01 sound power spectra repeatability. 
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Discrete point averaging consists of dividing your control surface 
into numerous discrete areas and measuring the sound intensity at the 
center of each discrete area (point) and summing over all areas. This 
technique offers maximum repeatability but at reduced spatially- 
averaged accuracy. The two techniques where both utilized initially 
and compared. The repeatability of the total radiated power using the 
discrete point technique, with 45 points, was 0.3 - 0.5 dB between 
successive runs. These repeatability values were selected as the 
criteria for the continuous spatial averaging technique. The results 
showed that the two techniques gave similar sound power spectrums and 
similar total sound power levels within an accuracy of 0.3 dB. For 
this reason, the continuous spatial averaging technique was selected 


due to its increased accuracy and reduced sampling time. 


The control surface was constructed of 1/4" balsa wood enclosing a 
total surface area of 0.542 square meters. Intensity spectra over the 
frequency range of 52-1600 Hz were taken and recorded on a B&K 2032 
Spectrum Analyzer. Microphone calibration was performed using a pis- 
tonphone and daily temperature and barometer measurements were taken 
to calculate the air density. To ensure that the two microphone were 
indeed phase matched, their phase was measured in a random noise field 
and compared. The microphones were placed at the open end of a long 
cylinder, equal distances from a speaker placed at the opposite end of 


the cylinder. This apparatus was similar in design to a calibration 





pistonphone, type 4200. Figure 2.4 shows the phase difference measured 
over a frequency range of 2560 Hz and clearly demonstrates their com- 


patibility over the frequency band of interest, 52 - 1600 Hz. 
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Figure 2.4: Phase difference between two microphones and their pre-amplifiers 
(channels А and B) 


In Section 5.2, the errors in dB resulting from these phase nis- 
matches are calculated and shown to be on the same order of magnitude 


as the inherent statistical fluctuations. 
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2.5.3 Power spectrum contamination 


Data below 52 Hz were ignored due to contamination by the 1,0,0 
room mode (38 Hz) of the semi-anechoic chamber. Also at this low 
frequency, the effectiveness of the anechoic treatment was reduced. 
The peak at 462 Hz and its harmonic at 924 Hz which is evident in all 
the spectra was due to the resonance of the cross-duct mode of the 
test section and not at all associated with the orifice source. For- 
tunately, the sound level and narrow bandwidth associated with this 
modal resonance contributes insignificantly to the total sound power 
caused by the orifice source. Frequencies above 1600 Hz contributed 
nothing to the total sound power and the spectral shape followed the 


same trend as from 1200 - 1600 Hz. 
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3. THEORY 


3.1 Introduction 


The development of a theoretical model to determine the scaling 
laws for orifice/louver flow noise was first accomplished in 1972 by 
Ffowcs Williams [8]. He treated the model problem of flow noise gen- 
erated by a body of turbulence near an infinite, rigid, thin screen 
with circular orifices. His assumptions were: a) the size of the ori- 
fice is much smaller than the distance between them, b) variations in 
pressure within individual apertures are neglected. While these 
assumptions are rather restrictive in their application to many indus- 
trial louvers, the theory does provide a fundamental result: that the 
sound radiated by orifices and louvers can be directly related to the 


mass flow driven through the aperture by the aerodynamic field. 


This fundamental result was used by Nelson [9] to develop a model 
for flow noise generated over practical perforated screens with circu- 
lar orifices. Due to the normal component of the turbulent fluctua- 
tions above the louver, flow is driven through the orifice resulting 
in a fluctuating volume flow through the upper surface. This fluctu- 
ating volume flow constitutes an equivalent monopole source at the 
upper surface. Since the volume flow into the upper surface of the 
orifice creates an equivalent outflow through the lower surface, there 


is an equal and opposite equivalent monopole at the lower surface. 
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The net field which radiates out of the louver and into the closed 
acoustic space is thus the combination of the lower monopole source 
and the contribution from the upper monopole source which is dependent 
on the degree of acoustical transparency of the louver. If the louver 
is acoustically opaque, then only the lower monopole source will radi- 
ate into the acoustic space. If the louver is acoustically transpar- 
ent, then the upper monopole source will radiate into the space, hav- 
ing a cancelling effect on the field of the lower monopole source. 
This would result in a net radiated field which exhibits properties of 


a dipole source. 


Nelson's equivalent source model gave an excellent description of 
the flow noise produced by perforated facing materials used in prac- 
tice even though he assumed a rather simplistic form for the driving 
aerodynamic velocity. His assumptions were: a) the wavelength for the 
sound is much larger than the orifice radius, b) there are many aper- 
tures per acoustic wavelength and no interactions of the flow occurred 
between adjacent apertures which can modify the aperture’s conducti- 
vity. In the sections that follow, the theory developed by Nelson and 
Ffowcs Williams will be utilized to model the flow noise created by 
turbulent flow over orifices and louvers. Once the equivalent source 
model has been developed, it will be used to collapse the experimental 
sound power spectra thus providing a verification to the theory's val- 


idity. The coordinated system used for this analysis is shown in 





Figure 3.1. 


Ks, 


Figure 3.1: Coordinate system for analysis 


3.2 The Equivalent Source Model for Orifices and Louvers 


Ffowcs Williams [8] and Leppington and Levine [10] developed an 
expression for the reflection coefficient, R, for a plane wave of 
single frequency, w, incident at a polar angle, 9, on an infinite 
louver. For low frequencies and orifice spacings that are large com- 


pared to the orifice size, the reflection coefficient, R, is 


R = 1/(1+4aiN/kcos#) (3.1) 


where a is the orifice radius, N the number of orifices per 


unit area of louver and k is the acoustic wavenumber, k = w/Co. 


The field produced by the lower monopole at yy can be determined 
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using the principle of reciprocity. The field at y, due to a distant 
source at xy, is given by p,(l*R) where p, is the incident plane wave 


striking the lower surface of the louver. Thus one can write 


P(yy) - (1*R) QGu). expCiklyy- xy D/75nlyy- xl (3:525 


where Q is the strength of the distance monopole at xy. By 
reciprocity, the field at xy due to the monopole at y, is given by 
interchanging the source and receiver position vector in equation 3.2. 
Equation 3.2 can be generalized to related the spectral density of the 
acoustic pressure field, Sp Ok. 9) to the spectral density of the 


source So(y,,#) such that 


S5, o) - |1«R|2* Sof 0) /16к212 (3.3) 


where r? = |xy-yy|?. Figure 3.2 is a plot of |R|? against the ratio 


(4aN/kcos6). 
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Figure 3.2: The infinite screen reflection coefficient as a function of the ratio 
(4aN/kcos@). From (9, Nelson, P.A., 1982). 





From this plot, it can be seen that, 


|R|? - 1 for (4aN/kcosé) < 0.3. 


This occurs for open area ratios on the order of 0.01 (see refer- 

ence 11). For all the orifice and louvers tested in this study, the 
open area ratio is less than 0.005 so that ||? - l. To verify this 
result for the reflection coefficient, since the assumption of many 
apertures per acoustic wavelength is not valid for the single orifice, 
a separate analysis was performed. Here, in two dimensions, the plane 
wave of single frequency, o, is incident normal on an infinite, rigid 
surface with a single opening, of dimension, 2b. Again the principle 
of reciprocity is used to predict the radiated pressure from the 
single orifice. Using matched asymptotic expansions, for b/A «« 1, 
the lowest order solution for the reflection coefficient is again, 


|R|? = l. The details of this analysis are shown in Appendix 1. 


Thus, 


S59) = SQ rk 0) /ат2г? (3.4) 


The instantaneous source strength q(yky,t) can be written in terms 


of the weighted integral of the unsteady normal component of velocity 





= 831 


over the area of the aperture. To take advantage of the homogeneous 
and stationary properties of the velocity in the plane of the wall, we 


can write 


qt) - - iwpf Је (ур, с) К(ѕр-ур) дА (Ѕү) (3.5) 


where y, is the vector position of some reference point on the aper- 
ture surface and K(sy-yy) is the weighting function, i.e., the 
function which converts the integral over the aperture area, s, = 
S],S2, to an integral over all space, defined as 
K(sk-yy) = 1 for sy-yy inside of the aperture area, A 
- 0 for sy - yy oustide A 
Following Corcos [12] the spectral density of the source strength, 
SQ(Yk +0) can be defined by 
SQ(yx,o) = р2027 fSy, (zy - sk o)K(zy) KC) dA Cz) dA Cs) (3.6) 
o œ 
where zy is a dummy space vector. Define ey by 
Ek = Zk = Sk 
thus 
$0(Ук›®) = plo? f fS, Ce KG Kee) dA CS) dA Ce) 


о cd 


= pul $ Si (eo 9)8 Ce) dACe) (3.6b) 


where the function, #(ey), defined as the orifice response function, 





is determined wholly by the geometry of the aperture and can be evalu- 


ated as 


a 


ô (ek) = ЈК (5р) К(ѕу+ер) ад (ѕу) = 1 х area of overlap 
„ср 


= (L-|<3|)(b-|<«|) £or |e3| < L, Jey] <b (3.6с) 
= 0 elsewhere 


for the rectangular shaped orifice shown below 
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Finally, we assume the cross-spectral density, Syy(tx,0), to have 


the product form 


Syy (ey, w) = Sy (w)A(Kez) B(Ke3) eth 4 (3.7) 


where S,(w) is the spectral density of the normal velocity fluctua- 
tions, K dis the convective wavenumber given by K = w/U, where U, is 
the convection speed of the boundary layer eddies and A(Ke4) and 
B(Ke3) are the amplitude functions for the normalized longitudinal and 


lateral cross-spectra, respectively. This form is similar to the 
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expression given by Corcos [12] for the cross spectral density of the 
boundary layer pressure fluctuations. The longitudinal cross-spectra 
of the eddies as they convect in the yj direction also contains the 
phase term exp(1Ke1). A(Ke1) is plotted in Figure 3.3 as a function 
of its argument for various ratios of separation distance to momentum 
thickness (r1/6*) for wall pressure. Since the wall pressure is 
related to the volume integral of the normal component of the velo- 
city, wall pressure cross-spectra data can be used to predict the lon- 
gitudinal and the lateral cross-spectra of the normal component of 
particle velocity. Comparing the wall pressure correlation length 
scales from [12] to the velocity correlation length scales from [17], 


we find that both the longitudinal and lateral length scales agree to 


within an order of magnitude, thus supporting the previous statement. 





Figure 3.3: Amplitude function for longitudinal cross-spectral density of wall 
pressure. From [16, B8N report No. 4110366, 1966] 





For the various orifice and louver samples tested, at a maximum 
value for Ei equal to l, A(Ke4) varies from a maximum value of l со 
a minimum value of 0.6 for orifice 02 at U-20m/s and f-250 Hz (the 
frequency at which the spectral level is 10dB down from the peak 
level). Throughout most of the samples tested, A(Ke]) - 1 over the 
frequency band of interest. Thus the longitudinal correlation lengths 
were typically several times the length of the orifice aperture. Thus 
the assumption is made that no decay in the eddy structure occurs in 


the longitudinal (yj) direction. 


The normalized lateral cross-spectra in the y4 direction is given 
by the term B(Ke3) which is assumed unity for no separation and to 
decay very rapidly with increased separation. To justify this assump- 
tion, B(Ke3) is plotted in Figure 3.4 as a function of its argument 
for various ratios of r3/8*. From this figure, the typical correla- 
tion length, 2, is calculated to be on the order of 2 cm. The ratio 
of the lateral correlation length, 2, to the aperture length, L, is 
0.18 thus supporting the initial assumption of rapid lateral decay. 
The correlation lengths are determined at the point where the ampli- 
tude function equals 0.4 and at the frequency corresponding to the 


peak power spectral level. 








Figure 3.4: Amplitude function for lateral cross-spectral density of wall pressure 
From [16, BBN Report No. 4110366, 1966] 


Substituting equation 3.7 into equation 3.6b and evaluating the 
integral over the rectangular shaped orifice, assuming Sy(yx,w) to be 


uniform over the orifice area gives: 


SQ(y,e) = p2w* S, (yy, o) +L f+ (2/K2) (1-cos (Kb) ) (3.8) 


where L is the orifice length in the y3 direction. See Appendix (3) 
for the details of the integration for the rectangular shaped orifices 
and reference [9] for the results for circular shaped orifices. 
Equation 3.8 states that the monopole spectral density is directly 


related to the normal velocity spectral density. 





3.3 Scaling Laws for the Radiated Sound Power 


Having determined the form for the source strength, we can now 
develop the scaling laws for the radiated sound power from various 
rectangular shaped configurations. Assuming that each orifice radi- 
ates sound incoherently, the total source strength for N apertures 
is simply N times the individual source strength for each aperture 
defined by equation 3.8. Substituting equation 3.8 into equation 3.4 


gives: 
Np2L2+S 1 Т? о 
Sp(Xk,9) = №125, (ур, о) (1-соѕ (КЪ) ) 0. /2я<г (3.9) 
The assumption is made that the convection velocity, U,, is equal 
to one half the freestream velocity, U,. Define the turbulent 
intensity, I, and the Strouhal number, St(4) and St(b) by: 
2 2 
(I^)y, 7 l/Us * J S,(yy, o) do (3.10) 
St(2) = f£2/U,, St(b) - fb/U, (3.11) 


Finally, integrating equation 3.9 over the bandvidth Ao with the above 


expressions substituted, yields 


mE 5 
(p(x, Jay = Ne?L- (1-cos(Kb))St(4) (I2) Us /2£nr* (amie) 


чш 





37 


For spherical waves, the sound power, P, can be determined by: 





(P(t))aw = NeL (1-cos(Kb))St(2)(12)q,, UU (3.13) 
Equation 3.13 is the final result for the flow noise created by tur- 
bulent flow over orifices and louvers with rectangular shape as 
described in Section (4.2). The controlling parameters for scaling 
sound power are thus: 

(1) Strouhal number based on the lateral correlation length, 2 

and aperture transverse dimension, b 
(2) the normal component of the Turbulent Intensity 
(3) the fifth power of the freestream velocity 


(4) the shape factor (l-cos(Kb)) 


In Section (5.10), this model is used to collapse the experimental 
sound power data for orifices 02 and 03 and louver Ll. It is clear 
from the collapse that equation 3.13 does indeed identify the predon- 
inant controlling parameters for rectangular shaped orifice and louver 
flow noise. Note, for circular shaped apertures, the above scaling 
factors remain the same except the shape factor is replaced by the 
function, Hj(Kd) where Hj is a Struve function of the first kind and 
d is the aperture diameter. See reference 9 for the details in cal- 


culating equation 3.8 for a circular aperture. 
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4. EXPERIMENTAL APPARATUS 


4.1 Semi-anechoic Wind Tunnel 


The experiments were carried out using the low turbulence, sub- 
sonic wind tunnel of the Acoustics and Vibration Laboratory at MIT. 
The wind tunnel details are shown in Figure (4.1) with the basic con- 
struction characteristics described by Hanson [13] and modifications 
described by Shapiro [14]. The wind tunnel is of open circuit design 
with the test section enclosed in a blockhouse. The tunnel consists 
of a set of flow straighteners, screens, a settling chamber and a 20:1 
contraction leading to a square test section 38cm by 38cm. The block- 
house walls, floor and ceiling were treated with urethane foam to pro- 
vide a semi-anechoic chamber within which acoustic intensity measure- 
ments were taken. Figure (4.2) details a cut-away front and top view 
of the chamber and test section. The acoustic treatment consisted of 
a base of 10cm foam, then a set of 10cm foam squares spread randomly 
and finally a set of 2cm foam squares also spread randomly. The 
effectiveness of the coating was determined by [14] using a horn 
driver and measuring the sound field decay with distance. The ane- 
choic treatment was effective between 200 Hz and 2000 Hz and only par- 


tially effective at 100 Hz. 
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Figure 4.1: Wind Tunnel Facility 


4.2 Orifice/Louver Samples 


Nine separate orifice/louver samples were tested. All samples 


were made of plexiglass and screwed into the sample mounting window 


which was mounted flush into the side to the test section as shown in 


Figure (4.2). The separate samples tested were: 


(1) 0-1: a 4 1/2" (height) by 1 2/3" (width) rectangular ori- 
fice, 3/4" thick with straight sided walls 


(2) 0-2: the above made of 1/4" thick plexiglass 
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Figure 4.2: Experimental set-up within the wind tunnel, front and top views 
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(3) 0-3: a 4 7/16" by 27/64" rectangular orifice, 1/4" thick 
with straight sided walls 

(4) 0-2A: 0-2 with 40 degree slanted walls, slanted away from the 
flow to reduce the downstream edge impact angle 

(5) 0-3A: 0-3 with 40 degree slanted walls 

(6) L-l: 4 rectangular orifices, each 4 7/16" by 27/64", separ- 
ated by 1/4", 1/4" thick 

(7) L-2: 54 circular orifices, each 27/64" in diameter, separ- 
ated by 0.6328" between centers, 1/4" thick 

(8) L-1A: L-1 with 40 degree slanted walls 


(9) L-2A: L-2 with 40 degree slanted walls. 


The louver configurations were designed to maintain the same total 
open area as the single orifices, 0-1 and 0-2. The 40 degree slanted 
walls were slanted to provide a blunt trailing edge and subsequently a 


sharp leading edge. 


4.3 Test Section 


The test section was made from 3/4" finished plywood, 38cm by 38cm 
square in cross section and was located on the centerline of the cham- 
ber. The test section was sealed and the ends were further covered by 
2" of urethane foam. The background power levels, measured with the 


orifice covered, were sufficiently low compared to the power levels 
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with the orifice exposed (see Section 5.3) to ensure signal to noíse 
ratios of at least 10 dB over the frequency range of interest and at 


all operating speeds. 
4.4 Turbulent Intensity Measurements 
Turbulent Intensity is defined by the relationship: 


2 e O 2л ( (4.1) 


Turbulent fluctuations were measured using a Dantec Constant Tempera- 
ture Anemometer with a 5 micron diameter tungsten wire. The output of 
the anemometer, being the total voltage E, was inputted directly to 
the B&K 2032 Spectrum Analyzer. Here the AC-coupled autospectrum, 
Gp(w), was measured. The mean-squared fluctuating velocity autospec- 


trum, G,(w) was determined from: 
(о) = (|H]%) Gelo) (4.2) 


where H is the transfer function. See Appendix 2 for the details in 
determining H. Turbulent Intensity spectra were measured at all four 
flow speeds for orifice 0-2 and louvers L-1 and L-1A. The position of 


the U probe was: 





(1) 


(2) 


Md. 


two-tenths of an inch into the boundary layer directly above 
the orifice at both the upstream and downstream edge (for the 
louver, the streamwise measurements were taken in the center 


of the furthermost upstream and downstream orifice, and 


half-way into the throat of the orifice and at the same 


streamwise location as described above. 
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5. EXPERIMENTAL RESULTS 

5.1 Introduction 

Figure 5.1 shows the geometry of the problem. Air flows at free 
speed U, over an orifice or louver which has been inserted flush into 
a long flat walled duct. The orifice was either rectangular, with 
streamwise dimension b and length L or circular with diameter d. 
The acoustic radiated power from the orifice/louver was measured at 
free stream velocities of 20, 30, 40 and 50 meters per second. The 
acoustic power was determined by sweeping the intensity probe over a 
control surface (a five sided boxed frame surrounding the orifice) 
then multiplying the resulting spatially and temporally averaged 
intensity by the total surface area of the control surface. Five ori- 


fices and four louvers were tested as described in Section 4.2. 


viver б уси : Ü 


upstream 
toS € 






(5 <ided box frame, 
Zach side dimension 32.9cm 
Figure 5.1: Schematic of orifice and probe geometry. Air flows at freestream velo- 
city, U, over an aperture of streamwise dimension, b. Sound intensity 
is measured by probe (microphone spacing, Ar) while continuously sweep: 
ing control surface. 
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To determine the contribution made solely by the orifice/louver 
source, two additional experiments were conducted. One experiment 
consisted of measuring the acoustic power radiated by the duct when 
the orifice was covered. The second experiment consisted of measuring 
the power flow through the orifice, in the absence of free stream 
flow, due to an artificially created acoustic field in the duct. This 
field was created using a speaker in the duct and was equal in magni- 
tude to that caused by the turbulent boundary layer at U = 50 m/s. 
This experiment determined the "open-window" contribution to the 
overall sound power. Finally, directivity measurements were taken for 
orifice 0-2 and louver L-2A at U = 25 m/s and the respective polar 


plots produced. 


5.2 Reactivity of the Sound Field and Sound Intensity Error Limits 


The reactivity of the sound field in the blockhouse (the closed 
acoustic space) was measured for various orifice and louver samples 
(see Figures 5.2 and 5.3). The Table below lists the measured dif- 
ference between the total sound pressure level, Lo and the total sound 
intensity level, L;, for the orifice 02A and the louver L2A at a 


radial distance of 2 feet. 





dB re 1Е-12 W/m? 








| | | | 
| Sample | — Reactivity (Lp:Li) | Us (m/s) | 
| | | | 
| | | | 
| 024 | 0.8 | 30 | 
| | 1.0 | 50 | 
| | | | 
| {2А | 0.9 | 30 | 
| | 1.0 | 50 | 
| | | | 
Table 5.1: Measured Sound Field Reactivity 
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Figure 5.2: Sound field reactivity for orifice 02A. Shown are the sound inten- 
sity and sound pressure spectral levels. The difference between the 
two is defined as the sound field reactivity. 
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Figure 5.3: Louver L2A sound intensity spectral levels versus sound pressure 


spectral levels (L; vs. Ly). The difference between L; and L 
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defined as the sound field reactivity. 


From Figure 2.4, the measured phase-mismatch for the sound inten- 


sity probe is: 


Frequency (Hz) Ýa (degrees) 
63 0.8 
500 0.2 
1000 0.5 
2000 0.9 


The maximum phase-mismatch between channels A and B for the Ana- 


lyzer is 0.3 degrees. 


With a spacer value of 50mm and the above sound 
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field reactivity, the phase reactivity nomogram of reference 6 gives 
the following value for the maximum phase of the sound field, between 


the two microphone positions. 


Frequency (Hz) $4 (degrees) 
63 e 
500 16 
1000 25 
2000 60 


Therefore the limits of the phase-mismatch are: 


Frequency (Hz) $¢°¢,°¢, (degrees) 6;+¢,*¢, (degrees) 


63 0.9 3.1 
500 15.5 16.5 
1000 24.2 25.8 
2000 58.8 61.2 


The above total phase-mismatch values are the limits of error in 
the measured phase at the frequency listed. The corresponding limits 


in dB for the measured intensity from the nomogram are: 


Frequency (Hz) Error limits (dB) 
63 + 2.5 
500 + 0.5 
1000 + 0.25 
2000 + 0.125 


The corresponding statistical standard deviation, e, of the mea- 


sured spectral levels about the mean value is: 
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є ~ 4.34 / /BT, , (dB) (5.1) 


where 


is the filter bandwidth and T, the total time record length. 


Substituting the appropriate values for B(2 Hz) and T, (150 sec) 


gives e = 0.25 dB. Thus the limits of error in measuring the sound 


intensity are approximately equal to the inherent statistical fluctua- 


tions (for frequencies above 100 Hz). 


5.3 Length and Frequency Scaling 


The three significant length scales for this experiment were: 


(1) 
(2) 


(3) 


the streamwise orifice dimension, b 
the turbulent boundary layer thickness, ô, above the orifice, 
and 


the lateral correlation length, £. 


Table 5.2 provides the applicable values for b, 6 and £ used in 


this experiment. 


The significant frequency scales were the non-dimensional strouhal 


numbers St(b) and St(42) defined by: 





$t(b/2) - f.b/U, (5:25 


St(2)  - o2/Uc - arf l/l 





| Orifíce/| | | | | $b | 6 | 
| Lower | beem) | Ш(т/$) | I ccm) |stem) | for 0-2 | foro-3 | 
| NE EC eL. ав 
| | | | | | | | 
| 0-4 | 4.23 | 20 | 2.4 | 4.19 | 0.99 | 3.95 | 
| | | | | | | | 
| 0-2,/ | 4.23 | 30 | 1.8 | 3.86 | 0.91 | 3.64 | 
| 0% | | | | | | | 
| | | | | | | | 
| 0-3,/ | 1.07 | 40 | 0.9 | 3.65 | 0.86 | 3.4 | 
| 034 | | | | | | | 
| | | | | | | | 
| L-1,/ | 1.07 | 50 | 1.4 | 3.49 | 0.83 | 3.29 | 
| 1L-1 | | | | | | | 
| | | | | | | | 
| L-2,/ | 1.07 | 50 | 2.1 | 3.49 | | | 
| v2 | | | | | | | 
| | | | | | | | 
| Full | 2.50 | 16.5| | 4.70 | 1.88| | 
| Scale | | | | | | | 
| | | | | | 


Note: Full scale assumes a louver positioned 15ft from the leading edge of a 
flat surface in seawater at 58 F. The louver consists of multiple rectan- 
gular shaped openings, each of transverse dimension, b=2.5cm. From Figure 
5.60 at é/b = 1.88, the Strouhal # = 0.19. Therefore the full scale 


frequency, fe, is: 
(с = (0.19)(16.5)/ (0.025) = 125 Hz 


Table 5.2 Experimental Length Scales 


where f. is the center frequency of the 1/3 octave band with the peak 
sound power level. Figure 5.4 is a graph of f, versus U, for orifice 


0-1 at a constant b. The data points reflect several separate 
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experiments measuring the sound pressure level at radial distances of 
9, 18, 24 and 36 inches from the orifice. Figure 5.5 plots f, versus 
U» for orifices 0-1, 0-2, 0-3 and louver L-l and L-2. Clearly there 

exists a linear relationship (or nearly linear) between f, and U, for 


a particular transverse dimension b. 
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Figure 5.4: Frequency/velocity relationship for orifice 01. Shown are four sep- 
arate experiments for various radial distances, r. 
Table 5.3 lists the non-dimensional Strouhal numbers for the above 
orifices/louver as well as the non-dimensional length scale, 6/b for 


U - 20n/s. The strouhal number, St-fb/U, has been defined to equal 


= 


E 


С 


м 





Center frequency of peak 1/3 OB (Hz) 


E 


that recommended by Demetz and Farabee [5] to compare results. Figure 
5.6 is the graph of the data in Table 5.3 along with the results of 


Demetz and Farabee (5] for TBL flow over orifices enclosed by solid 


cavities. 
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Figure 5.5: Frequency/velocity relationships for orifices 01, 02, 03 and louvers 
11,12 
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Strouhal number. 
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| | | | 
| Orifice/ | St(b) = fb/U | Non- dimensional | 
| Louver | Strouhal $ | length scale 3/b | 
| | | | 
| | | | 
| 01 | 0.17 | 1.22 | 
| 0-2 | 0.23 | 1.22 | 
| 03 | 0.12 | 4.79 | 
IB 1 | 0.076 | 4.79 | 
| L-2 | 0.048 | 4.79 | 
| | | | 


Table 5.3: Won-dimensional Frequency and Length Scales at U = 20 ms. 
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Figure 5.6: Strouhal number versus boundary layer/aperture width ratio 
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The decreasing trend in the Strouhal number with increasing 6/b 
for the two cases agrees reasonably well. Having a knowledge of the 
boundary layer thickness, one can enter Figure 5.6 at the appropriate 
value for &/Ъ ага determine the Strouhal number for orifice or louver 
flow noise. From the Strouhal number and the geometry of the apera- 


ture(s), one can predict f.. 


5.4 Structural Radiation Due to TBL Excitation 


To determine the contribution from structural radiation caused by 
the excitation from the TBL, an experiment was conducted at all free 
stream velocities with the orifice/louver covered. Identical sound 
power spectrums were recorded as described in Section 4.1. The 
resulting background power spectrums levels, Figures 5.7, 5.8, 5.9 and 
5.10, were 10 to 46 dB lower than with the orifice exposed. From 
these figures, the major contribution from the orifice/louver sources 
was in the range of 65 to 1000 Hz. Background noise was within 15 dB 


of the orifice/louver noise below 50 Hz above 1000 Hz. 
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Figure 5.7 Signal to noise for orifice 01 at U = 30m/s. 
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Figure 5.8: Signal to noise for orifice O2 at U = 40 ms. 





OB RE IE-12 WATT 


POWER. 


DB RE 1E-12 WATT 


POWER, 


58. 


40. 


38. 


cp. 


18. 


= 


18. 


ШҮ 


p. SBB. 1882. 1580. 
FREQUENCY (HZ) 
Figure 5.9: Signal to noise for orifice 03 at U = 20 avs. 
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Figure 5.10 Signal to noise for louver Li at U = 50 ws. 
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5.5 "Open Window" Contribution 


Previous work by Spinka [15] determíned the acoustic pressure 
field within the test section at various free stream velocities. To 
determine the contribution from the open window affect, a similar 
acoustic pressure field was generated within the duct using a speaker 
driven a B&K random noise generator, type 1402, and a B&K band pass 
filter set, type 1612. The resulting power which radiated through the 
orifice was measured as described in section 2.5 and compared to the 
power radiated by the orifice source. Table 5.4 lists the resulting 
power levels for both cases at U - 50m/s. The conclusion was that the 
open window affect results in levels which were typically 23dB lower 


than the levels due to the orifice source. 


| | | | | 
| 1/3 08 center | Open windon | Orifice 0-1 | | 
| frequency | power in 1/3 08 (2) | power in 1/3 08 (3) | (3-2) | 
| | | | | 
| | | | | 
| 400 | 53.8 | 77.6 | 38 | 
| 500 | 48.2 | 72.7 | 24.5 | 
| 630 | 45.9 | 69.2 [M2 | 
| 800 | 43.7 | 67.9 | 2.2 | 
| | 1000 | 46.2 | 66.4 | 20.2 | 
| 1250 | 41.0 | 64.5 | 3.5 | 
| | | | | 


Table 5.4: "Open window" power levels versus orifice O -1 power levels at U = 50m/s 
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5.6 Orifice Source Power Spectra and Total Power 


The measured sound power spectra versus the free stream velocity, 
U-20,30,40 and 50 m/s, for 0-1, 0-2, 0-3, L-1 and L-2 are shown in 
Figures 5-11, 5-11a, 5-1lb, 5-11с, 5-114. The total sound power radi- 


ated by each orifice/louver source at the above speeds is given ín 


Table 5.5. 
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Figure 5.11 Orifice 01 sound power spectral level versus free stream velocity 
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Figure 5.11(a): Orifice 02 sound power spectral levels versus freestream velocity 
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Figure 5.11(b): Orifice 03 sound power spectral levels versus freestream velocity 
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Figure 5.11(c): Louver L1 sound power spectra levels versus freestream velocity 
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Figure 5.11(d): Louver L2 sound power spectral levels versus freestream velocity 





| | | | 
| Orifice/ | Free Stream | Total Sound Power | 
| Louver | Velocity U(m/s) | (dB re 1 E-12 watts) | 
| | | | 
| | | | 
| 0-1 | 20 | 63.3 | 
| | 30 | 73.3 | 
| | 40 | 79.6 | 
| | 50 | 84.0 | 
| | | | 
| 0-2 | 20 | 62.7 | 
| | 30 | 73.2 | 
| | 40 | 78.4 | 
| | 50 | 83.0 | 
| | | | 
| 0-3 | 20 | 62.9 | 
| | 30 | 72.2 | 
| | 40 | 76.3 | 
| | 50 | 79.7 | 
| | | | 
| L-1 | 20 | 74.5 | 
| | 30 | 82.5 | 
| | 40 | 89.2 | 
| | 50 | 94.5 | 
| | | | 
| L-2 | 20 | 64.9 | 
| | 30 | 73.6 | 
| | 40 | 80.1 | 
| | 50 | 85.3 | 
| | | | 
| 0-2А | 20 | 62.0 | 
| | 30 | 71.8 | 
| | 40 | 77.2 | 
| | 50 | 82.2 | 
| | | | 
| L-1A | 20 | 58.1 | 
| | 30 | 68.4 | 
| | &0 | 73.3 | 
| | 50 | 7?.7 | 
| | | | 
| L-2A | 20 | 55.3 | 
| | 30 | 65.0 | 
| | 40 | 69.7 | 
| | 50 | 76.3 | 
| | = | 


Table 5.5: Orifice/louver total sound power levels versus freestream velocity,U, 





Several observations are made from reviewing the above data. 


5.6.1 The effect of louvering on sound power 


(a) louver, L-l, radiated the most acoustic power, averaging 10dB 
more at all speeds than the single orifice though the total 
open areas were the same (see Figures 5.12 and 5.12a). 

(b) louver, Ll, showed strong resonant peaks and had the narrow- 
est half-power bandwidth of all the samples tested. 

(c) louver, Ll, radiated significantly more power than louver L2 
in the lower frequency bands but approximately equal power in 


the middle to high frequency bands (see Figures 5.13 and 


5.14) 
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Figure 5.12 The effect of louvering on power, Li versus 02 at U-30n/s 
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Figure 5.128: The effect of louvering on sound power, Li versus 02 at U=50m/s 
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Figure 5.13: Louver geometry effect on sound power, L1 versus at L2 at U = 30m/s 
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.14: Louver geometry effect on sound power, Li versus L2 at U = 50m/s 


louver, L2, had the broadest spectra of all the samples 
tested with equal open area. 

louver, L2's total radiated power was within 2 dB of orifice 
02's power in all speed ranges but the spectral shape was 
significantly broader. As shown in Figures 5.15 and 5.16, 
louver L2 radiated significantly more power than 02 in the 
middle to higher frequency bands (10 dB more on average) but 
slightly less (approximately 1-2 dB) power in the lower 


frequency bands. 
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Figure 5.15: Effect of louvering on sound power, L2 versus 02 at U = 30m/s 
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Figure 5.16: Effect of louvering on sound power, L2 versus 02 at U = 50 m/s 
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Thus the effect of louvering, simply by dividing the total open area 
into more numerous, smaller areas does not by itself reduce the total 
radiated flow noise and will most probably result in increased radi- 


ated sound power in the higher frequency bands. 


5.6.2 Orifice throat thickness effect on radiated power 

Figure 5.17 shows the orifice throat depth effect on sound power 
by comparing orifice 01 to orifice 02. Orifice 02 radiated more power 
(approximately 3dB) in the higher frequency bands but maintained 
approximately 1-2 dB less power in the more dominant lower frequency 


bands. 
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Figure 5.17: Orifice depth effect on sound power, 01 versus 02 at U = 30 and 50 m/s 
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Thus the throat thickness was more predominant in effecting the 
spectral shape than reducing the total radiated power. The effect of 
reducing the orifice throat thickness appears to broaden the spectral 


shape while maintaining the same total radiated sound power. 


5.6.3 Streamwise dimension, b, effect on sound power 


The spectral shape of orifice 03 had the widest half-power band- 
width. The effect of reducing the streamwise dimension, b, was to 
broaden considerably the power spectra while reducing only moderately 
the total sound power. From section 3, the controlling shape factor 
was (l-cos(Kb)). Thus the predicted result from reducing b by four 
fold would be a subsequent increase by four fold in the power spectra 
bandwidth. Reviewing again the power spectra for orifices 02 and 03 
(Figures 5.11a and 5.11b), we see that indeed the power bandwidth for 
03 is approximately four times that of 02 (defining the bandwidth to 


be between the one-quarter peak spectral power points). 


5.7 Speed/Power Laws 


The speed/power law is of the form: 


P= 6." 


where C is a value which is independent of speed. The data in Table 
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5.5 was utilized to determine the total radiated power dependence on 
speed. The equations used to determine the value of the speed expo- 


nent, n, were: 


P4/P5 - (Uz/U2)” 
10 log (P3/Pre¢) - 10 log(Po/Pr.r) = n 10 log(U3/U5) 


Lp,” Lp," n log(Us/U2) 
n7 (Ly,- Lp}/(10 log (U3/U2)) 


where the subscripts 3 and 2 correspond to 30 and 20 m/s. Similar 
equations hold for U,/U5, Us/U2, U¿/Uz, Us/Uz and Us/U¿. Table 5.6 
lists the calculated value for n for each sample tested. The values 
for n varied from 3.3 to 6.8 with a mean value of 4.9 and a standard 
deviation of 0.3. Excluding the extremely low values of n recorded 
for 03 at the higher free stream velocities, the mean value of n was 


5.0 with a standard deviation of 0.2. 


The speed/power relationship for orifices 01, 02 and louvers Ll 
and L2 is plotted in Figure 5.18. Also shown is the least squares fit 
to the data points satisfying the equation: 


(Lp, - LP4) = n 10 1og(U,/U4) (5.2) 


From the fit, the best value for n is 4.97. 
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Figure 5.18: Speed/Power relationship for four samples (PzCU") with least squares fit 
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Thus the concluding dependence is: 


Total radiated power increases with the 5th power of the free 


stream velocity 


5.8 Directivity of orifice Ol and louver 12А 


Directivity measurements from orifice 01 апа louver L2A were 
taken at U - 25m/s. Figure 5.19 details the coordinate system used 
for the measurements. The measurements indicated the source was 
nearly omni-directional. Power levels in the XZ (azimuth) plane dev- 
lated by 0.9 dB from the maximum level. In the YZ (polar) plane, power 
levels deviated by 2.6 dB from the maximum level. Off the major axis, 
power levels deviate by 3.6 dB from the maximum level directly above 
the orifice source. Table 5.7 lists the measured sound intensity in 
dB re 1E-12 watt/m? at a radial distance of r - 27 inches and the 
various polar and azimuth angles. Figure 5.20 is the azimuth direc- 
tivity plot for orifice 01 and louver L2A and clearly supports the 


conclusion of a nearly omni-directional source. 
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Figure 5.19 Coordinate system for directivity measurements. Reference is the 
aperture center. Polar angle, 6, measured from 2 axis, azimuthal 
angle, ¢, measured from -x axis. Measurements described by spheri- 
cal coordinates, r,0,4. 


| | 
| Total Intensity | 
| Polar angle Azimuthal angle dB re 1E-12watt/me | 
| | 
| Orifice 0-1 | 
| | 
| 0 90 57.3 | 
| 0 45 56.9 | 
| 0 15 56.4 | 
| 0 135 57.1 | 
| 0 165 56.8 | 
| 0 105 57.2 | 
| 0 75 57.0 | 
| -85 90 54.9 | 
| -45 90 56.5 | 
| 0 90 57.5 | 
| 45 90 57.0 | 
| 85 90 55.9 | 
| -45 30 54.7 | 
| -45 25 53.9 | 
| -45 35 55.0 | 
| | 
| | 





Table 5.7 Directivity measurement for orifice 0-1 and louver Lea 
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Table 5.7 Directivity measurement for orifice 0-1 and louver L2A (Continued) 


Total Intensity 
dB re 1E - 12watt/n* 
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5.9 Edge Angle Effect on Radíated Sound Power 


5.9.1 Introduction 


Earlier it was shown that indeed the orifice source provided both 
the spectral character and the predominant amount of the total power. 
The strength of this source is due in part to the eddy interaction 
with the downstream edge of the orifice. To determine the relative 
importance of this edge interaction compared to the other significant 
parameters identified in Section 3, three additional orifice/louver 
samples were tested. They were orifice 0-2A and louvers L-1A and 
L-2A. The details of these samples are described in Section 4. The 
only difference between these samples and the previously tested 
samples was the leading and trailing edges were slanted away from the 
flow by 40 degrees. This provided both a blunt trailing edge where 
the eddies strike and a sharpened leading edge (essentially reducing 


the upstream aperture to 1/32", the radius of the edge). 


5.8.2 Angled edge orifice/louver power spectra 


The measured sound power spectra for the above samples are shown 


in Figures 5.21, 5.22 and 5.23. The total power radiated by each 


sample at the various freestream velocities is given in Table 5.5. 
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Figure 5.21: Orifice 02A power spectra level versus freestream velocity 
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Figure 5.22: Louver LÍA power spectral levels versus freestream velocity. 
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Figure 5.23: Louver L2A power spectral levels versus freestream velocity. 


Several observations are worth noting: 


(a) the effect of angling the edges is more predominant for the 


(b) 


louvers than the single orifice. 


the effect on louver L-1A is a 14-16 dB reduction in total 
radiated sound power. The large resonant peaks which occurred 
in the spectra for L-l are not present in the spectra for 
L-1A. Figures 5.24 and 5.25 compare the two louvers at U=20 
and 50 m/s. Not only are the resonant peaks eliminated by 
angling the edges but also the radiated power is reduced by 


5-7 dB in the higher frequency bands. 





DB RE IE-12WATT 


POWER, 


POWER, OB RE IE-12WATT 


60. 


32. 


40. 


38. 


20. 





502. 1202. 1588. 2090. 
FREQUENCY (HZ) 


Figure 5.24: Edge angle effect on sound power spectral level, L1 versus LÍA 
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Figure 5.25: Edge angle effect on power spectral level, L1 versus L1A at 


U = 50 m/s 





(c) Figure 5.26 compares the power spectral levels of L-2 and 
L-2A. The result for this louver is a 8.6 to 10.4 dB reduc- 
tion in sound power equally divided across the frequency 
bands. 

(d) for the single orifice, O-2A, the effect of angling the edges 
resulted in only a 0.7 to 1.4 dB reduction in sound power. 
Figures 5.27 and 5.28 compare the power spectra for orifice 
0-2 апа O-2A at U-20 and 50 m/s. Here again the effect of 
angling the edges resulted in a reduction in sound power in 
the higher frequency bands but in this case, essentially no 


change in the lower frequency bands. 
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Figure 5.26: Edge angle effect on sound power spectral level, L2 versus L2A at 
U = 50 m/s. 
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Figure 5.27: Edge angle effect on sound power spectral level, 
02 versus 02A at Û = 50 m/s. 
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Figure 5.28: Edge angle effect on sound power spectral level, 


02 versus 02A at U = 20 m/s. 
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Combined effect of louvering and angling edges 


Taken together, the effect of louvering and angling edges results 


in a reduction of 3.4 to 7.5 dB in radiated sound compared to the 


single angled orifice (both of equal total open area). This reduction 


most probably is caused by two factors: 


(a) 


(b) 


the louvers shed eddies which after convecting downstream to 
the trailing edge are the same order in size as the edge 
thickness. This results in a more beneficial interaction when 
the eddies impact the angled edge. The orifice, however, with 
a four times larger aperture, sheds eddies which are much 
larger in size upon impact with the trailing edge and are 


therefore less affected by the thin angled edge. 


louvering, by reducing the transverse dimension, b, interferes 
with the longitudinal correlation of the shed eddies. This 
interference will be more predominant in the lower frequencies 
where the longitudinal correlation length is larger than the 
aperture transverse dimension, b. As shown in Figures 5.29, 
5.30, 5.31 and 5.32, the reduction in sound power is predom- 
inantly in the lower frequency bands. There is, however, more 
radiated power in the angled louver ín the hígher bands by 


approximately 1-3 dB. 
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Figure 5.29: Effect of Louvering on sound power spectra level, L2A 
versus 02A at U = 30 m/s 
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Figure 5.30: Effect of louvering on sound power spectral level, LIA 
versus 02A at U = 30 m/s 
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Effect of louvering on sound power spectral level, 12А versus 02A 
at U = 50 ms 
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Effect of louvering on sound power spectral level, LIA 
versus 02A at U = 50 m/s. 





Thus one can conclude that effective louvering does reduce the 
total radiated power by reducing single mode resonance and by inter- 
fering with the lower frequency, longitudinal eddy correlation. This 


results in effectively broadening the power spectrum. 


5.9.4 Louver geometry effect on sound power 


As with the straight walled louvers, the multiple circular aper- 
ture louver was quieter than the multiple rectangular aperture louver 
by 1.4 to 3.6 dB in total radiated power. As shown in Figures 5.33 
and 5.34, louver L-2A had approximately the same spectral shape as 
louver L-1A but with considerably less power in the low frequency 
bands. Again this reduction can be attributed to the additional 
interference of the lateral correlation due to the reduced dimension 
in the aperture’s height. To test this conclusion, one can calculate 
the frequency where the lateral correlation length, 2, equals the 
aperture diameter, d for louver L-2A and at freestream velocity of 50 
m/s. From Figure 3.4, at B(K(r4))-0.4, the value of Kr4 equals 1.2. 
Setting r4 equal to 0.01 m (27/64") and U-50 m/s gives a limiting 
frequency of 450 Hz. Thus one would predict diminishing reductions in 
sound power for frequencies above 450 Hz. Reference to Figure 5.34 
shows that, indeed, at frequencies above 450 Hz the effect of lateral 


correlation interference on reducing sound pover is diminishing. 
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Figure 5.33: Louver geometry effect on sound power spectral level, 
LÍA versus L2A at U = 30 m/s. 
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Figure 5.34: Louver geometry effect on sound power spectral level, 
LIA versus L2A at U = 50 m/s. 
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Figure 5.35: Turbulent intensity spectral levels for louver L1 versus freestream 
velocity. 


5.10 Describing the Velocity Field Surrounding the Samples 


In Section 3, the radiated sound power was shown to be directly 
related to the normal component of the turbulent intensity. Several 
measurements of the turbulent intensity spectrum were taken as 
described in Section 4.4. The assumption was made that the turbulent 
intensity in the streamwise and normal component directions were 
equal. Thus 3 dB was subtracted from the measured intensity spectral 
levels to give only the normal component of the turbulent fluctua- 


tions. The following observations are made: 





2 бт 


(a) the turbulent intensity measured 2 tenths of an inch into the 
boundary layer was relatively unaffected by the probe loca- 
tion in the transverse direction and the various orifice 
geometries. Figure 5.35 is the measured normal component of 
the turbulent intensity for louver L-1, with the probe at the 
center of the upstream aperture. These spectral levels were 
used in the following section to collapse the sound power 


spectra. 


(b) there were significant reductions in turbulent intensity 
(approximately 3 dB ) between louvers L-1 and L-1A when meas- 
ured at the leading edge of the leading aperture. This reduc- 
tion in intensity supports the measured reduction in total 
sound power and also supports the conclusions that the sound 


power is related to the driving aerodynamic field. 


Despite a fairly rigorous attempt to maintain the block-house air 
tight, a finite and measurable pressure differential existed between 
the blockhouse and the duct static pressure. This pressure differen- 
tial resulted in a small amplitude, mean flow entering the ori- 
fice/louver from the closed acoustic space. The magnitude of this 
transverse mean flow was greatest in the louver samples and measured 
on average, fourteen percent of the freestream velocity. The exact 


contribution of this mean flow to the radiated sound power in not 


БЕЙТ © 


known but it 1s estimated to be of minimal impact since the turbulent 
Intensity of the entering flow measured several orders of magnitude 


less than the normal component of the freestream turbulence. 


5.11 Verifying the Scaling Laws 


In Section 3, an equivalent source model was developed for rectan- 
gular shaped aperture orifices and louvers. To ascertain the validity 
of this model, the sound power spectra for orifices 02, 03 and the 
louver Ll were collapsed by normalizing the spectral levels by the 
magnitude of the equivalent source model. The resulting collapsed 
Spectra are shown in Figures 5-36, 5-37 and 5-38. It is evident by 
the collapse that the equivalent source model adequately describes the 
controlling parameters for orifice and louver flow noise. The predic- 
tion of the actual magnitude of the sound power is fairly accurate for 
the samples with the smaller apertures but is too high by about a fac- 
tor of 100 for the larger aperture openings (01 and O2). The probable 
reasons for this are (1) overestimating the effective aperture width, 
b, (2) the estimation for the lateral correlation lengths and (2) the 
assumption of equal components of turbulent intensity. As mentioned 
in Section 3, for orifice 0-2, the longitudinal correlation amplitude 
function, A(K(e4) does attenuate to a value of approximately 0.6 for 
the slower speeds indicating that some eddy decay occurred in the 


aperture. For this case, one could assume a smaller transverse 





length, b', over which the eddy decay was minimal and use this length 
scale for evaluating the source strength. Defining b' at the poínt 
where the magnitude of the amplitude function equaled 0.95 results in 
an effective aperture length of 0.5 inches versus the actual dimension 
of b-1.67 inches (for orifice 0-2 at U-20m/s and f-250 Hz). This 
results in a reduction of 4.2 dB ín the predicted spectral level and 
could account for some of the error in predicting the larger aperture 


radiated sound power. 
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Figure 5.36: Scaling law verification for orifice 02 





18. 
CD 
D 
E ð. 
> 
м 
"~ 
al 
ww 
4 
o -1B, 
Pr 
"~ 
D 
x 
euh. 
0 
О 
d -2@. 
G E 
3 m 
X Y) 
CQ М 
= -30 ۲ 
a 
us 
m 
Ó 
J 
© 
-AQ. 


1QLOG (P- FC/NpL (1-Coe (Kb) ) St (1) IU), dB 


о. вов 9.858 


8. 108 8.158 В. 2080 8.250 9.3900 


STROUHAL NUMBER BASED ON b, Stsfb/U 


Figure 5.37: 
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Scaling law verification for orifice 03 
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Scaling law verification for louver L1 
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6. CONCLUSIONS 


The acoustic source created by turbulent flow over orifices and 
louvers was investigated by measuring the radiated sound intensity, 
utilizing an intensity probe, at various wind speeds in a semi- 
anechoic wind tunnel. From the intensity spectra, the radiated power 
and the directivity of the source was determined. The goal of this 
research was to describe the source and identify the controlling par- 
ameters which affect its spectral shape, radiated power and directi- 
vity. The orifice and louvers investigated have geometries similar to 
those used in industry in both air and water environments. The 
effects of various leading/trailing edge geometries were also investi- 


gated. 


The source of noise has been modelled on the basis of the theoret- 
ical analysis by Ffowcs Williams, Nelson and Corcos. The equivalent 
source of noise was modelled as two equal but opposite monopole 
sources located in the upper and lower aperture of each orifice. The 
strength of each source is determined by the mass flow driven through 
the apertures by the normal component of the turbulent fluctuations. 
The resulting radiated sound field, both into a closed acoustic space 
and radiating away from the orifice, is dependent on the acoustical 


properties of the louver. 
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To ascertain the validity of the source model, two orifice and one 


louver power spectra were collapsed by normalizing the spectral levels 


by the magnitude of the equivalent source model. It is evident by the 


collapse that the equivalent source model adequately describes the 


controlling parameters for orifice and louver flow noise. The con- 


trolling parameters for the geometries tested were: 


(1) 


(2) 


(3) 
(4) 


made: 


(a) 


(b) 


a 


Strouhal number based on the lateral correlation length, 2, 
and aperture transverse dimension, b 

the normal component of the Turbulent Intensity, 

the fifth power of the freestream velocity and 

the shape factor (1-cos(Kb)) where K is the convective 


wavenumber, K = w/Uc, and Uc is the eddy convection velocity. 


result of this research, several important conclusions are 


sound intensity measurements, using an intensity probe, can 
successfully describe an acoustic source, even though 
measurements are taken within the source's near field. 

the errors associated with measuring intensity, using an 
intensity probe, can be as small as the inherent statistical 


deviation of the random fluctuations. 





Ge) 


(d) 


(e) 


(£) 


(g) 


EO. 


the equivalent source model developed to describe the behav- 
ior of orifice and louver flow noise successfully describes 
the controlling parameters and provides valuable insight for 
predicting the actual magnitudes of the radiated sound spec- 
trum. 

the 1/3 octave band containing the most acoustic power from 
orifice and louver flow noise can be predicted from the non- 
dimensional Strouhal number based on the aperture's trans- 
verse dimension and the non-dimensional length scale, 5/b 
where § is the turbulent boundary layer thickness. 

for aerodynamic flows, the orifice/louver source provides 
both the dominant amount of total sound power and the spec- 
tral character of the radiated noise. Contributions due to 
structural radiation from the ducting and ducting “open win- 
dow" effects are not significant. 

the effect of louvering, simply by dividing the total open 
area into more numerous, smaller areas does not by itself 
reduce the total radiated flow noise and will most probably 
result in increased radiated sound power in the higher 
frequency bands. 

the effect of reducing the aperture throat thickness is to 
broaden the spectral shape of the radiated power spectrum 
while maintaining approximately the same total radiated sound 


power. 





(h) 


(1) 


(j) 


(k) 


(1) 
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the effect of reducing the aperture's transverse dimension 
can be adequately described by the controlling shape factor 
in the equivalent source model, (l-cos(Kb)). 

directivity measurements show a nearly omni-directional 
source in the on-axis planes with off-axis power levels 
reduced by 2-3 dB. 

taken together, the effect of louvering and angling the lead- 
ing and trailing edges resulted in reductions to sound power 
of 3.4 -7.5 dB compared to the equivalent area, single aper- 
ture, angled orifice. 

the measured reductions ín sound power due to effective lou- 
vering can be explained by (1) reductions in the normal com- 
ponent of the turbulent intensity due to the eddies impacting 
a blunt traíling edge and (2) destructive interference with 
the eddy longitudinal and lateral correlations. 

the louvers with numerous circular apertures and angled edges 
radiated the least acoustic power with a very broad, flat 


power spectrum. 
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APPENDIX 1: MATCHED ASYMPTOTIC EXPANSION SOLUTION TO A RADIATING 
ORIFICE IN AN INFINITE RIGID WALL 
The two dimensional problem to be solved is the far field radiated 
pressure at xy, created by a driving aerodynamic field at yy, within a 
finite size aperture in an infinite, rigid wall. Using the principle 
of reciprocity, the above field can be determined by first solving for 
the sound field at yy created by a distant monopole at xy. Then, the 
desired solution is obtained by simply replacing the position vectors 
of the far field monopole and the aperture source. The solution to 
this problem will be analyzed using the principles of Matched Asymp- 


totic Expansions (MAE). 


The coordinate system used in the analysis is shown below, 


a 


% 


with boundary conditions, 3p/dya - 0 for |y,| = b and Y = O. 


The governing equations for a plane wave of single frequency, o, 


at normal incidence are: 





С) А 


у2р + kp - 0 (A1.1) 
DS NOE) 
p = Ae (Al.2) 
where k is the acoustic wavenumber, k = w/cg. Also define the 


non-dimensional terms, f - ot, yj - Ку], yo = куд, p= D/A апа є = kb 


for the parameters in the outer problem. 


Substituting the above outer parameters into equations Al.1 and 
Al.2 gives the governing equations and boundary condition for the 


outer problem 
у2р + р - 0 (А1.3) 


i(y2-7) 


p-e (A1.4) 


To lowest order, «20, and the orifice closes. Thus the outer 


problem becomes 


P-61p] * ó62po * *** as e-0 (A1.5) 


with 6) = 1. From the boundary condition, we can determine pj to be 


ly? -1y9 
р = руі + рү = е Е ‚ уз < 0 (А1. 6) 





egg 


suppressing the e" Í7 factor, where Pi». Pr are the incidence and 
reflected waves, respectively. Note that the boundary condition 1s 


satisfied 


др\ ly? -iyo 

iem ce -e ) = 0 for yo = 0 (Al.7) 

ду2 

Choose also the condition that pz = O for y2 > 0. To open the 
orifice, choose the inner variables y1 - yı/e and y2 = yo/e. 


In inner variables, 


p = exp(i(eyo-r)] + exp[-i(eyy-r)] 


exp(-1r)(exp(ey2) * ехр(-єуд)) 


2l dy _ a 
- exp(-ir)(leeyo - ——2 - ++++1- уо - 2 +666) 
2 2 
_ ey : 
- 2e7 (1- ——? « Ket)) as e 9 0,y2 fixed (yo - 0) (A1.8) 
2 


The governing equation with the inner parameters substituted into 


equations Al.3 becomes 
v2p + «2р -0 (A1.9) 


with boundary conditions 3p/9y2 = 0 for yy = O and |yı| > 1. Let 


p be defined by 


р ~ До(є)ру + Ду (є)рә + eee (A1.10) 





EON 


for the inner expansion. The resulting inner problem coordinate sys- 


tem is 





Since V2p4 = 0, solve the inner problem by conformal mapping, Pi = 
A Rew. Using a solution of the form, iz - cosh w, (suppressing the 
eit factor) we can determine 
1 13 
ы ~ In 2 z + i«/2 + - tees as |z|-.o (A1.11) 


2 
2х22 2x4x4z 








where z = rei?. 


Since 2 = z/e, 


w - fn 2z - fne + in/2 + @(e2) for e > 0, z fixed (zoe) (А1.12) 


Therefore, utilizing equation A.11, we can determine Ag(€) Pj to be 


ло(є)ру ~ Abol fn 2r - ne + in/2 + O (e2)) (A1.13) 
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Now matching the leading terms in the outer and inner solutions, 
equations Al.8 and Al1.13, leads со До = + 1//п є апа А = - 2. Then 
ö62(e) = 1/!n e in the outer expansion. p» satisfies the Helmholtz 
equation with k = 1 and must behave like ln r forr + 0, « + 0 andr 


fixed. Clearly 


(1) 
ро = € Hg (т) (A1.14) 


satisfies this requirement and also satisfies the rigid wall condi- 
tion. As r э 0, 
(1) 21 21 


Hg (r) ~ — inr * 1 * — (4 - 1fn2)++.. 
x x 


where y is defined as Euler's constant. Thus to match with the inner 
solution, A- - 2- C « 2i/r, C must be equal to ir. Therefore the 
final outer expansion becomes 

Y2  -Y2  ixH(D(r) 


P- e11 (e + е + — ot oee ) (A1.15) 
fne 


where the first term in brackets is the incident wave, the second term 
is the reflected wave and the lowest order orifice effect is the last 
term which describes a monopole type radiation. Substituting for the 
non-dimensional terms and applying the principle of reciprocity, one 


obtains the radiated sound from a 2 dimensional, single slot, line 
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source 
~ kx2 -kxo irH( (x) -iot 
P- Ale + е pr o e (A1 16) 


Ankb 


= De — i(r-n/4) 
As r +o, н M + EE е So the for field radiated 


sound for the 2D, single slot, line source becomes 


_ кх) -kx2 X TEENS C 
P - Ale te * Ж е Je (A1.17) 


pte 


Equation A1.17 provides the result that the far field radiated 


pressure for a 2D, single slot source to lowest order is composed of: 
(1) the incident and rigid wall reflection 


(2) 2D line source having a l/r dependence (cylindrical spread- 


ing). 


Thus one can infer that it is plausible for the 3 dimensional, 
single orifice source to have a reflection coefficient, |R|? = 1 for 


values of b/A << 1. 
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APPENDIX 2: MEASURING TURBULENT INTENSITY 
In this appendix, the transfer function to relate the mean and 


fluctuating component of the hot wire anemometer output voltage to the 


mean and fluctuating component of the particle velocity is developed. 


The basic relationship between output voltage, E, and particle 


velocity, U, is defined by King's Law 
E = A + BUP (A2.1) 
Now define the instantaneous voltage and particle velocity by 


E = Eg * e', le'| «« |Eol (A2.2) 


U = Ug + u’, lu’ | << 00| (А2.3) 


where Eg is the mean of E, Ug is the mean of U, e' and u' are the 
fluctuating components whose expectations are zero. Substituting 
equations A2.2 and A2.3 into equation A2.1, expanding and dropping the 


higher order terms of u' and e' gives 


ES + 2Ege' = А + BU (1+ [nu' /Ug]) (42247 


Taking the expectation of equation A2.4 gives 





OS 
E? - A 4 Bun (A2.5) 
0 0 


Subtracting equation A2.5 from equation A2.4 gives the desired 
relationship between the fluctuating components of the anemometer 


output voltage, e’, and the particle velocity, u’, 


l-n 
u' - ([2bEgUg ]/Bn)e' (A2.6) 
Thus, defining the term in parentheses as H, the resulting rela- 
tionship between the spectral densities of the fluctuating velocity, 


Su(f), and the output voltage, Sp(f), is 


Su(£) = |u|? sgc£) (A2.7) 
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APPENDIX 3: MATHEMATICAL STEPS IN DETERMINING Sq Gio 9) 


From equation 3.7, the cross spectral density is assumed to be of 


the form 


Suv Cio 9) = S,(w)AlKeı)B(Kez3) © exp (ike) (A3.]) 


As described in Section 3.2, the normalized longitudinal cross- 
spectra amplitude function, A(Kej), is identically equal to one for 
all samples except Ol and 02 at the lower freestream velocities, where 
at a minimum it equals 0.6. So assuming no decay in the eddy struc- 
ture in the yj direction (i.e., A(Ke1) equals one) and S,(yy,o) to be 
uniform over the aperture surface, the source strength spectral den- 


sity becomes 


So(yk,9), - ^u S Or 9) f BO) exp Ie) 8 Ce) dA C) (A3.2) 


where «1, єз аге the position coordinates as shown below. 
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Substituting for 6(ey) we have 


ike 
1 
SQ(yk,e) = pw Sy(yk.») Î f B(Ke3) e  (L-|<3|)(b-|«1|) (43.24) 
со o 
* dA (ey). 
Since the lateral correlation й, is typically one-tenth or less 


the length of the aperture, the integral in the ez direction can be 


evaluated, independent of the actual limits of the integrand, as 


| B(Ke3) (L-|e3])de4 = L2 - #037) 


(A3.3) 
со 
where 
| B(Keg3)deg = £ (АЗ . За) 
a 
Evaluating the integral over the e] directions gives: 
iKe, 2 
e (b-[eı |)dej - =z (1-cos(Kb)) (A3.4) 
K 


со 
Therefore the final result, neglecting the 1? term is: 


2 
SQ(Yk+0) = SOLL) (1-cos(Kb)) (A3.5) 
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The remaining term needed to fully describe the magnitude of the 
source strength spectral density is 2. The lateral correlation length, 
2, is determined at the point where the amplitude function, 

B(Ke3) equals 0.4 and the frequency equals the center frequency, f,, 


as defined in Section 5.3. 
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